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Abstract. The problem of galaxy formation and its dependence on thermodynamic properties is addressed by us- 
ing Eulerian hydrodynamic numerical simulations of large scale structure formation. Global galaxy properties are 
explored in simulations including gravitation, shock heating and cooling processes, and following self-consistently 
the chemical evolution of a primordial composition hydrogen-helium plasma without assuming coUisional ion- 
ization equilibrium. The galaxy formation model is mainly based on the identification of converging dense cold 
gas regions. We show that the evolution at low redshift of the observed cosmic star formation rate density is 
reproduced, and that the galaxy-like object mass function is dominated by low-mass objects. The galaxy mass 
functions are well described by a two power-law Schechter function whose parameters are in good agreement 
with observational fits of the galaxy luminosity function. The high-mass end of the galaxy mass function includes 
objects formed at early epochs and residing in high-mass dark matter halos whereas the low-mass end includes 
galaxies formed at later epochs and active in their "stellar" mass formation. Finally, the infiuence of two other 
physical processes, photoionization and non-equipartition processes between electrons, ions and neutrals of the 
cosmological plasma is discussed and the modifications on galaxy formation are examined. 
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• I— I 1. Introduction 

X 

^ Processes like gravitation, shock heating, radiative cool- 
ed , ing, photoionization and non-equipartition, among others, 
play a crucial role in the evolution of the thermodynamic 
properties of baryonic matter. As galaxies originate in cold 
and dense gas regions, any change in the gas thermody- 
namic properties should have an effect on their formation. 
This connection is the issue here addressed. Numerical 
simulations have the significant advantage of being able 
to include a large set of physical processes involved in 
galaxy formation. Even if a phenomcnological descrip- 
tion of this process needs to be adopted, global galaxy 
properties are now extensively studied, like the cosmic 
star formation rate density, the galaxy mass/luminosity 
function, the clustering properties, etc. Such results use 
semi-analytical approaches (,So mervi ll e fc Pr imack 1 999 ; 
ICole et al hybrid approach es, combination of semi- 

analytical and N-body methods i Kauffmann et al]ll999t 
iHat ton et al. 2003}; and hydrodynamical N-body simula- 
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tions using Smooth-Hydrod ynamic-Parti cle JPearce et all 
l200lt 'Weinberg et al.'2002^. Lagrangian (|GnedirJll996() or 

Eulerian (Gen 1992) methods. These complementary com- 
putations are all based, but now in a sophisticated way, on 
the fundamental ideas that galaxy formation results from 
the g as accretion and i ts coo l ing (.Sunya e v fc Zcrdovicl^ 
197i iRees fc OstrikMl Il977t ISilkl il977t Iwhite fc ReeJ 
Il97q) into a network of overdense structures created by 
the gravitating dark matter on the large scales of the uni- 
verse. Moreover, quantities computed from numerical sim- 
ulations can now be compared with the huge amount of 
avail able observational da ta. The galaxy luminosity func- 
tion ( Binggeli et alJll98'8l) is widely estimated, in different 
surveys, towards fainter magnitude, in several wavelength 
bands, and for different classes of galaxies (Kladgwicki 
2002 ). Neve r theles s some discrepancies remain: Fig. 1 in 
iCross et al.l l|200ll) shows a dispersion of a factor of 2 at 
the characteristic luminosity, L*, and a factor of 10 at 
O.OIL,. A number of studies attempt to retrieve the galaxy 
mass function from the galaxy luminosity function by us- 
ing stellar population synthesis models. Numerical Simula- 
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tions take the opposite approach: their first output is mass 
and by using stellar population synthesis models they can 
provide luminosity. Here we focus, among other things, 
on the cosmological evolution of the galaxy mass function 
and examine in detail its shape at 2: = and z = I. We 
then derive mass-to-luminosity ratios to compare with the 
observed galaxy luminosity function. 

As this paper is the second in a series of three, 
focusing on the influence of thermodynamics on galaxy 
formation, the purpose is not to compute sophisticated 
models of galaxy formation but rather emphasize the 
description of cosmological gas. Hence to keep the model 
free of parameters as much as possible, we only consider 
in the simulations the dominant physical processes: grav- 
itation, shock heating, radiative cooling, but neglect e.g. 
feedback processes. A model of galaxy formation is also 
introduced. The first part of the paper examines global 
galaxy properties, the cosmic star formation rate density, 
the galaxy mass function and the epoch of formation. 
The main results are the following: 1. Galaxy formation 
is a hierarchical process mainly driven by the amount of 
available cold gas in the inter-galactic medium; 2. The 
majority of the high-mass galaxies form at early epochs; 3. 
The galaxy population at any given redshift is dominated 
by a significant fraction of low-mass galaxies formed at 
early as well as late epochs. The halo dark matter mass 
function is also explored and a preliminary study of the 
galaxy distribution inside halos is presented. The second 
part of the paper analyzes separately the influence of 
photoionization, from ultraviolet background radiation, 
and the influence of non-equipartition processes between 
ions, neutrals and electrons of the cos mological plasma. 
Non-e quipartition has been scrutinized in ICourtv fc Alimil 
l)2004|) (Paper I) using two numerical simulations: the 
first one taking into account non-equipartition processes 
and denoted by S^t and the second one, denoted by 
Sn, in which equipartition between species is forced. 
The former simulation allows each species to carry its 
own internal energy whereas the latter one assumes that 
ions, neutrals and electrons have the same temperature. 
That paper concludes that a significant fraction of the 
inter-galactic medium (the plasma inside gravitationally 
bound structures), that is accreted in not too dense 
structures and at temperatures in the range 10^-10^ K, 
is out of equilibrium and warmer in the S^t than in the 
5'iT simulation. Non-equipartition processes are likely to 
be dominant before the end of the reionization epoch. As 
galaxies are accreting their gas from the inter-galactic 
medium in the temperature range 10'*-10^ K, this implies 
an infiuence of the non-equilibrium thermodynamics on 
the galaxy formation process. Quantifying this change 
is one of the purposes of this paper. The third paper in 
this series (Courty & Alimi, in preparation), will quantify 
how galaxy clustering properties and the cosmological 
bias are modified. 

This paper is organized as follows. Numerical simu- 
lations and the galaxy formation model are described in 



section|51 Section|31presents galaxy properties: the cosmic 
star formation rate density, the galaxy-like object mass 
function, and the epoch of formation. The dark matter 
halo mass function and the galaxy distribution inside the 
biggest mass halos are discussed. Fits of the galaxy-like 
object and dark matter halo mass functions are given. 
The influence of photoionization and non-equipartition 
processes are shown in Section 0] and Section |S1 respec- 
tively. Conclusions are given in Sectional 



2. Numerical simulations 

The simulations were performed with a 3 dimensional 
N-body/hydrodynamical code, coupling a Particle-Mesh 
method for computing gravitationa l forces with a 
Eulerian method (Tevssic r et alJll998|) . The simulations 
include shock heating, radiative cooling, photoionization 
processes, non-equipartition processes between the ions, 
electrons and neutrals of the cosmological plasma, and 
galaxy formation. The features of the simulations, ana- 
lyzed here, are the following: the GO and Gl simulations 
include shock heating, radiative cooling and galaxy 
formation. They only differ in the resolution. In addition 
to these processes, the GP simulation includes photoion- 
ization processes and GNE includes non-equipartition 
processes but not photoionization. We refer to paper I for 
details about shock heating treatment, non-equipartition 
processes and the radiative cooling terms. These latter 
terms include collisional excitation, coUisional ionization, 
recombination ^, bremsstrahlung and Compton scatter- 
ing. To use the same notations as in paper I, GO, Gl and 
GP are 5'iT simulations, in which equipartition between 
species (ions, neutrals and electrons) is forced and the 
cosmological plasma has a single temperature. The GNE 
simulation is a Syr simulation with each species having 
its own internal energy. 

The GP simulation takes into account the ionization 
and heat input from an ultraviolet background radiation 
to reproduce conditions after the reionization epoch. The 
photoionization and heating rates are computed from the 
evolution of the hydrogen and helium densities and from 
the spectrum of the ultraviolet background radiation J{v). 
The radiation is considered a spatially uniform field over 
the computational volume. The density evolution equa- 



^ The helium recombination rate of these simulations has 
an incorrect temperature dependence, although the expres- 
sion for oljj++ in Table 1 of Paper I is commonly used in 
the literature. The correct expression should involve the nu- 
clear charge of the helium atoms (seeiSpitzcr : a„++ = 

3.36 X 10""' (1 + (Te6/4)" '^)"^re"^/^r^3°-2. Since collisional ex- 
citation is the dominant net cooling term this mistake should 
have a limited effect on the results. 
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tions, Eq. (5) to (7) in Paper I, are now: 



dt 
dt 
dt 



where n-no, nH+, n-HoO> "-Hc+i "-Hc++ ^-^^ "-e are the six 
densities of the primordial composition hydrogen-helium 
cosmological plasma, and Pi and a; are the ionization and 
recombination rates (given in Table 1 in Paper I) . Photo- 
ionization rates Pi are expressed by: 



hv 



ai{v)dv 



with i denoting the species H", He°, He+, Ui is the effecti ve 
cross-section for species i (taken from lOsterbrockl ljl989l) ). 
hvi is the ionization energy for the species i, and J{v) is 
the background radiation intensity. The photoionization 
processes are also a heating source (see Eq. (4) in Paper 
I) and the heating rates are expressed by: 



H, 



AttJ{v) 
hv 



{v){hv — hvi)dv 



(2) 



The shape of the background intensity spectrum is 
defined by a function F{z) characterizing the evolu- 
tion with redshift of the ultraviolet background radiation 
l)Katz et al.lll996|) : 



J{u)=F{z) (- 



-1 



erg cm ^ s ^ sr ^ Hz ^ 



(3) 



where Hvh — 13.598 eV is the hydrogen ionization thresh- 
old. Since this evolution is little-known at high red- 
shift, we estimate J(i^) from observational measurements 
llScott et al.l I2n02^ and from numerical works fjGnedh] 
l200n|) . This includes the decline of the ultraviolet back- 
ground radiation intensity observed between z ~ 1 and 
z = 0, the sharp increase before z ~ 6 and the shallow 
evolution before z ~ 7. We start the reionization at red- 
shift 10.5 with the bulk of the transition between z = 7 
and z — 6. The function F{z) is plotted in Fig. ^ and its 
expression, in units of Jq = 10" 
is: 



F(z) 



' lQ(Allog(l+z)+Bl) 
'^Q(A2log{l + z)+B2) 
lQ(A3log{l + z)+B3) 

Jo(4/l-t-z)4-" 
Jo 

[Jo(2/(l + z))-3 



11.5 > l-fz > g 
8 > 1 + z > 7.7 
7.7> 1 + z > 7 

7 > 1 + z > 4 
4 > 1 + z > 2 
2 > 1 + z 



(4) 



with Al = -7.27,^2 = -98.32, v43 ^ -29.52, Bl = 
-19.28, B2 = 62.94, B'S = 1.94. 



Fig. 1. Evolution with redshift of the ultraviolet back- 
(1) ground radiation intensity. 



As it is useful for the discussion in this paper, we now 
compare the gas distribution in the three kinds of simu- 
lations. Figure |2] illustrates the baryonic mass fraction in 
temperature-density diagrams; the top and bottom pan- 
els are extracted from Paper I. The top panel is computed 
with a simulation including hydrodynamical shocks and 
radiative cooling. The middle panel includes photoioniza- 
tion processes and shows that they are dominant over 
cooling processes only in low and middle dense regions 
(see also IWeinberg et all lll997j) '): low density regions are 
heated up to temperatures between 10'^ and less than 10^ 
K. The competition between photoionization heating and 
cooling due to adiabatic expansion results in the concen- 
tration of the gas on a slope T = Tn ( p/pV with Tn ^ 6.10^ 
K and 7 around 0.6 iKatz et all ()l996V 'H ui fc Cxuedh] 
()l997t)). The bottom panel displays the gas distribution in 
the GNE simulation and shows a non- negligible warm gas 
fraction. As the influence of non-equipartition processes is 
larger at high redshift, the isocontours are plotted at z = 5 
(see discussion in Paper I) . 

Although differences exist in the low-density regions, 
altogether the three gas distribution diagrams show a 
common feature: a peak in the high density, cold region. 
This gas being the reservoir for galaxy formation, we thus 
present in the first part of this paper galaxy properties 
in a simulation only including the dominant processes 
involved in galaxy formation, namely hydrodynamical 
shocks and radiative cooling. The influence of photoion- 
ization and non-equipartition processes will be considered 
in separate sections. 

We now turn to the description of the galaxy forma- 
tion model. Numerical simulations of large scale structure 
formation currently do not allow for the formation of ob- 
jects beyond the scale of a grid cell, a few times 10^ M0, 
this mass being much larger than the mass of a single 
star. The problem is bypassed by considering the physical 
conditions needed to form a galaxy. The most important 
condition is that the gas cloud is collapsing, meaning that 
the cooling time is less than the dynamical time or the 
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GP z = 



lO" r 



10= r 



^ 10^: 



10^ 
10^ 




0.1 1. 10. 10^ lO-" 10* 

Pb/ <pb> 



GNE z=5 




1.0 10.0 100.0 

Pb/<Pb> 

Fig. 2. Isocontours of the baryonic mass fraction per in- 
terval of the baryonic density contrast and per interval of 
temperature computed for the Gl and GP simulations at 
z = and for GNE at z = 5 (top, middle and bottom 
panels, respectively). The increase in the mass fraction 
scales from dark to light. Note that Gl and GNE are, 
for these plots, computed without galaxy formation. Solid 
curves on the top panel show the ratios tcooi/tff = 0.1, 1, 
10. 



free fall time ijRees fc Ostriked HqTtI) . The top panel of 
Fig. |2] illustrates in a temperature-baryonic density dia- 
gram that dense and cold gas regions are located inside 
the iso-contour tcooi/^ff = 1- The galaxy formation model 



then consists of the identification of the gas satisfying this 
criteria. A fraction of the baryonic matter is turned into 
a "stellar" particle describing the amount of stellar mass 
produced during the process of galaxy formation. Galaxy- 
like objects are then defined by a collection of this "stellar" 
m ass (the term "galaxy -like object" has already been used 
in lEvrard et alJ l|l994|) . although describing high-density 
contrast baryonic clumps). To make sure that gas regions 
giving birth to galaxies are correctly identified we add 
other criteria described below. The purpose of this paper 
being to show how modifications of the gas thermody- 
namics have an effect on the galaxy formation process, we 
then deliberately keep the number of free parameters low. 
This model, although very simple, gives consistent results 
between the properties of the galaxy-like objects in the 
simulation and their observational counterparts. 

To express the condition tcooi < iff, we define the 
cooling timescale ^cooi, computed from the internal energy 
variation of the gas E/E, and the dynamical time or free 
fall time: 



te 



32Gp 



(5) 



One of the other conditions requires that the size of the 
gas cloud must be less than the Jean's length given by: 



1/2 



(6) 



Note that the total matter density, including dark mat- 
ter, baryonic matter and "stellar" particles, is used in the 
expressions of the dynamical time and the Jean's length. 
It is clear that the Jeans criterion is reliable only down 
to the mass resolution, since the size of gas clouds itself 
is limited by the spatial resolution of the simulation. A 
third condition is that the gas must be in a converging 
flow: V ■ V < 0. Finally the baryonic density contrast, 
Sb = {Sp/p)b, must be higher than a threshold (1 -I- Sb)s- 
This is taken to be the value of the baryonic density con- 
trast at the turnaround, 5.5, compute d in th e top-hat col- 
lapse spherical model ( Padmanabhanl ll99,'^ . 

To estimate the amount of "stellar" mass formed, we 
express the variation of the baryonic mass as the ratio be- 
tween the available baryonic mass tob and a characteristic 
timescale i,: 

dt ~ U 

The integration of this expression on a timestep At ^ t—to 
gives: 



niBit) = mB(io)exp(- — ) 



(8) 



where TOB(io) is the baryonic mass initially present. The 
"stellar" mass formed is then: 



= mB(to) - "^B(^) — nT-B(io) 



At 



(9) 



Then in each cell, checking the four criteria described 
above, a fraction of the gas is turned into a "stellar" 
particle. Each of these particles carries its mass m* and 
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Simulation Lbox dx Mdm 

(h-^Mpc) (ft-^kpc) (Mq) 



Mbm 

(Mo) 



GO 32. 125. 2.01 x 10* 3.09 x 10^ 

Gl, GNE 16. 62.5 2.51 x lO'^ 3.87 x 10^ 

GP 11. 43. 8.17 x lO** 1.25 x lO** 

Table 1. Parameters of the simulations. Lbox is the co- 
moving length of the computational volume, dx is the spa- 
tial resolution of the grid, Mdm is the mass of the dark 
matter particle and Mbm is the initial baryonic mass en- 
closed within a grid cell. 



its epoch of formation given by the scale factor a,. The 
mass m* is computed using Eq. Q with me (to) the bary- 
onic mass enclosed within the grid cell at each timestep 
and with the characteristic time — max(iff,10^ yr). 
Figure lA.ll in Appendix fXl shows the evolution in rcdshift 
of the "stellar" particle mass function for the different 
simulations: the "stellar" particle mass ranges between a 
few times 10^ and ~ 2.10® Mq. The "stellar" particles are 
involved in the computation of the gravitational poten- 
tial and their evolution is treated in the same way as the 
coUisionless dark matter. 

At any redshift two catalogs of objects are created: 
one consisting of dark matter halos and one of galaxy-like 
objects. Halos and galaxies are defined by grouping 
either dark matter particles or "stellar" particles with a 
Friend-of-Friend algorithm. This algorithm joins together 
all particles separated by a distance proportional to 
the link parameter rj. We take rj = 0.2. We exclude 
from the dark matter halo catalog groups with less than 
10 particles. This threshold is denoted Afmin- But the 
galaxy-like object catalog is allowed to include objects 
with a lower Mmin, meaning that each "stellar" particle 
is considered a galaxy-like object. The influence of these 
two parameters on the galaxy mass function, 77 and Mmin, 
are discussed in Appendix IbI 

Unless otherwise stated, the results of this pa- 
per are given for a A— cold dark matter model (A — 
CDM). The parameters of the simulations arc: Hq = 
70. km s-i Mpc"\ = 0., f^m = 0.3, f^A = 0.7, 
l^b = 0.02/1-2 with Ho/ 100. The initial density fluc- 
tuation spectr u m use s the transfer functions taken from 
Bardccn ct al. ("1986") with a shape parameter given by 
Sugivama ( 1995 . ). The fluctuation spe ctrum is normalized 



to COBE data ijBunn fc Whitelll997t) leading to a filtered 
dispersion at i? = 8 h-^ Mpc of as = 0.91. The number 
of dark matter particles is Np — 256"^ and the number of 
grid cells is iVg = 256'^. Three computational volumes are 
used, described, as well as the simulation parameters, in 
Table [H 
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Fig. 3. Evolution with redshift of the cosmic "stellar" 
mass formation rate density for the Gl (dashed line), 
GO (dot-dashed line) and GP (dotted line) simulations. 
The observational data of the star formatio n rate density 
are ov erplotted for o ur cosmology: Ha data:lGallego et"al] 
1I1995I) (filled circle).lTresse fc Maddo^ ll 19981) (filled trian- 
gle). 'Yan et"ai] l)l999j) (fille d star), [H opkins et al.1 l|200(t 
(hollow diamond); UV data : IConnoUv et al. (1997) (open 
square ). fTrever et al.l (Il998l) (open tri a ngle) .ISteidel et alJ 
11999) (op en circle). ISullivan e t al.' ("20001) (diamond); 
FIR data: iRowan- Robinson et al.. (,1997.) (heavy cross) , 
iFlores et all (ll999^ (cross); 1.4 GHz d ata: ICondonI (Il989t) 
(filled sq uare). 'Haars ma et al.l ()2000|) (six-pointed star), 
Serjeant et al. (2002 ) (open star). UV data and data by 
Gallego et air()l995|) are corrected for extinction. 



K o 
^ o 



^ 1 



20 



l+z 



Fig. 4. Same as Fig. Obut for a SCDM scenario, for the 
Gl (dashed line) and GO (dot-dashed line) simulations. 
The observation points are the same as described in Fig.|21 
but overplotted for a SCDM cosmology. 



3. Galaxy properties 

3.1. Cosmic "stellar" mass formation rate density 

Figure O shows the redshift evolution of the "star" forma- 
tion rate (SFR) density for the GO and Gl simulations. We 
compute the amount of "stellar" mass formed per year and 
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per unit of volume. The SFR density shows a strong de- 
crease at low redshift and peaks around redshift ~ 3 in GO 
and 2.5 in Gl. However, the ratio between the SFR at the 
peak and the SFR at z = 6 is lower than between the SFR 
at the peak and at z = 0, SFR{z = Zpcak) /SFR{z = 6) ^ 
3.2 against SFR{z = z^^^],) / S F R{z = 0) - 10 in the GO 
simulation (these values become 2.8 and 6.3 in Gl, respec- 
tively). This trend is related to the hierarchical nature of 
cold dark matter models. At high redshift baryonic mat- 
ter is accreted in dark matter potential wells and cools to 
form galaxy material. Then low-mass structures merge to- 
gether to form larger mass units. Additional gas accreted 
is then shock heated towards higher temperatures leaving 
insufficient time for the gas to cool and condense. This 
results in the decrease in the star formation rate at low z. 

Another illustration of the influence of the large scale 
structures on galaxy formation comes from the compar- 
ison with a different cosmological scenario. We run the 
same simulations but now computed for a standard cold 
dark matter model {SCDM) with the following param- 
eters: Ho = 50. km s-i Mpc"\ Qk = 0., flm = 1., 
Qa = 0., f^b = 0.02/1-2 ^i^j^ ^ ^ i/o/100. As in the 
A — CDM scenario, the fluctuation spectrum is normal- 
ized to COBE data, giving cg — 1.1. For computational 
reasons the SCDM Gl simulation is only performed up 
to z = 1 and the GO up to z = 2. The density fluctua- 
tion power spectrum has more power on small scales than 
the A — CDM model, resulting in a larger amount of cold 
gas at high redshift. Differences between these two mod- 
els are strikingly illustrated by the cosmic star formation 
rate density (Fig. 0J. The higher amount of dark matter 
in the SCDM , resulting in a slightly larger erg than in 
the A — CDM scenario, creates deep potential wells at 
high redshift and as a result the peak of the SFR density 
is reached at z = 6, or even higher z, leading to a steep 
decrease at low redshift. 

Figure 121 shows that the SFR density depends on the 
spatial resolution of the simulations: the highest resolu- 
tion simulation, Gl, has a SFR density amplitude higher 
at low redshift than the GO simulation. Decreasing the 
box length allows to include in the computational vol- 
ume density fluctuations with lower wavelengths, result- 
ing in a higher mass fractio n of available co l d gas (sim- 
ilar trends are discus sed in IWeinberg et alJ l)l999|) and 
lAscasibar et al.l l)2002 ^ ) . Figures and |B| display, for the 
same simulations as GO and Gl but without galaxy for- 
mation, the evolution with redshift of the baryonic mass 
fraction in different ranges of temperature, corresponding 
to the main phases of the inter-galactic medium: the "dif- 
fuse" phase with T < 9.10"^ K, the "cold" phase with tem- 
perature in the range 9.10'^-2.10^ K, the "warm" phase in 
the range 2.10^-5.10^ K and the "hot" phase at T > 5.10^ 
K. The increase in the amount of cold gas with decreas- 
ing redshift is less dramatic for the GO simulation: between 
z = 5 and z = 0, the cold gas mass fraction goes from 5.2% 
to 26% whereas it goes from 4.6% to 41% in Gl. However 
changing the resolution does not change the evolution of 
the SFR density and we are less interested in making quan- 




Fig. 5. Evolution with redshift of the baryonic mass frac- 
tions computed in different temperature ranges, for the 
GO simulation but without galaxy formation: T < 9.10^ 
K ("diffuse", dot-dot-dashed line), 9.10^ < T < 2.10-* K 
("cold", dashed line), 2.^ < T < 5.10^ K ("warm", dot- 
dashed line), T > 5.10^ K ("hot", dotted line). 




Fig. 6. Same as Fig.|31 but for the Gl simulation without 
galaxy formation. 



titative estimates than analyzing qualitatively the process 
of galaxy formation. 

Observational data of the star formation rate density 
are overplotted in Fig. |3| These data show a great dis- 
parity and the amplitude of the SFR density is expected 
to change due to dif ferent observational bias corrections 
I Hopkins et al.l200lh . Nonetheless a strong decrease is ob- 
served at low redshift. Moreover observations do not cur- 
rently agree whether the SFR density peaks after z = 1 or 
reaches a plateau at higher redshift. Both star formation 
rate densities in our simulations reproduce the slope at 
low redshift and the amplitude obtained in the Gl sim- 
ulation is consistent with observations (contrary to the 
SCDM scenario showing a decrease at earlier epochs). 
This agreement and the previous discussion point out that 
galaxy formation is mainly driven by the amount of avail- 
able cold gas enclosed in gravitationally bound structures, 
the strongest constraint of the galaxy formation model 
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lyiiiR in the condition tcooi < tg (see also IPearce et al.l 

3.2. Galaxy-like object mass function 

The rest of Section |3 is now devoted to the description 
of our two galaxy catalogs and examines an important 
property of galaxy formation, the galaxy-like object mass 
function. Discussions about their epoch of formation and 
their location inside dark matter halos will follow. 

Figure 13 shows the cosmological evolution of the mass 
function of the galaxy-like objects in the GO simulation. 
The catalog covers a wide range a mass, from 10^ M0 up 
to 10^^ Mq. As we do not know if the decrease in the range 
M < 10^ Mq is a numerical effect (see also Murali ct alj 
we discuss in the following the shape of the mass 
function for M > 10^ Mq. At any redshift the comoving 
number density of objects per bin of mass increases as 
the mass decreases. This trend is different from the stellar 
particle mass functions (Fig. lA.lT) . The galaxy mass func- 
tion shows a clear evolution from z = 5 to z = 0: more 
and more objects of high-mass are created, increasing the 
mass range of galaxies towards the high-mass end. Since 
galaxy formation results from a hierarchical process, the 
mass of bigger objects increases with decreasing redshift, 
shifting the knee of the mass function towards larger mass. 
On the other hand the low-mass end of the mass function 
shows an increase as redshift decreases to only z — 2. At 
lower redshift the trend is inverted with a decrease in the 
number density per bin of mass of low-mass objects. 

On the whole the galaxy mass function shows a 
characteristic shape: a strong decrease at the high-mass 
end beyond a characteristic mass, preceded by a shal- 
lower slope in the intermediate mass range, between 
10^-10^^ M0, and a steeper slope at the low-mass end up 
to 10® M0. Note the sharp transition at z = around a 
few times 10® Mq, shown more clearly in Fig.|HI 

Figure [7| presents one of the main results of the paper: 
the galaxy mass function is significantly dominated by a 
low-mass galaxy population, M < 10^'^ Mq, whatever the 
redshift is. This population can be linke d to the observe d 
faint luminos i ty ga laxy population l|Lovedavl l)l998l) : 
iNorberg et al.l l)2002l) . and references therein). In fact this 
faint population covers a wide variety of galaxies: galaxies 
with or without emission lines llZnr.ca et a] .111 9971^ . galax- 
ies chara cterized by a si gnificant star formation activity 
l)Lin et al. 1996: Lovedav et alJll999(l dwarf galaxies of 
morp hological and spectral late types llMarzke et alJl99'^ 



1^2^), low surface brightness galaxies fSuravberrv et al.| 
1997t) . blue compact objects ( Guzman et al, 1997). The 
discussion below about the epoch of formation of the 
galaxy-like objects will show that, in the simulations also, 
the low-mass end of the galaxy mass function includes 
an inhomogeneous population. This population of faint 
galaxies is important for galaxy evolution, since by their 
number they are likely to contribute to the cosmic star 
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Fig. 7. Galaxy-like object mass function at different red- 
shifts for the GO simulation. 



formation rate density. Their relationship to the environ- 
ment and their implication in galaxy mergers still needs 
to be addressed. 

The galaxy luminosity function is generally fitted by 
a standard Schechter function combining a power-law 
and a n exponential function at the bright end jSchechteil 

[Tizi): 



*w ^ £ U: 



exp - — 



(10) 



where a is the slope of the power-law, is the charac- 
teristic luminosity at the break, and 0* is a normalization 
parameter. TablejUlists parameter fits, extracted from the 
literature, of the galaxy luminosity function for different 
surveys. In order to differentiate mass from magnitude, the 
symbol M is used for magnitude in this table and through- 
out the paper. As pointed out in the introduction, discrep- 
ancies remain between surveys, enlightening the fact these 
surveys cover different galaxy populations. Observational 
results show that the galaxy luminosity function is differ- 
ent for galaxy populations selected by color, morphology 
or environnement. The 2dF survey, for instance, shows 
how the shape of the galaxy luminosity function changes 
when galaxies are divided according to their star forma- 
tion activity IjFolkes et al...l999il . 
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Fig. 8. Analytical fits (solid curve) of the galaxy-like object mass function (cross) for the GO and Gl simulations at 
z — and z = 1. Fits are computed using Eq. Hll|l . 



Most of these surveys use a single power-law Schechter 
function and the a slope is generally found to be in the 
range —0.9 and —1.2. A number of estimates of the galaxy 
luminosity function conclude that it cannot be described 
correctly, in the entire range of luminosity and specially 
at Mb > —15, by a single power-law Schechter func- 
tion l)Lovedavlll998|l . Two works, reported in Table [3 
use a non-standard Schechter function to parameterize 
the galaxy luminosity function: a two power-la,w Sche chter 
function described below, is used i n iLoved av' ( 1997) and 
two analytical functions are used in lZucca et al.i (|1997|), a 
standard Schechter function and a power-law at the faint 
luminosity end. Both report slopes in the faint luminosity 
range of —1.82 and —1.57, respectively. 

The mass function of the galaxy-like objects in Fig. 
clearly suggests similar conclusion. Therefore to account 
for the low-mass galaxy population we choose a two 
power-law Schechter function now combining a standard 
Sche chter furiction with a p power-law at the low-mass 
end l|Lovedavlll997l) : 

jm = m:[m:J + J ^^^^ 

where (t>{M) = diV/dAf is the numerical density of objects 
per interval of mass and per unit of volume, (/>*, M» and 



a are analogous parameters as used in Eq. l(TU|l . and Mt 
is the transitional mass between the two power-laws. The 
mass function of the galaxy-like objects is fitted, in the 
range M > 10® M©, using a least square method weighted 
by the mass function itself to ensure the statistical relia- 
bility of the result. Figure |S1 and Table |31 give fits and 
parameters determined at z = and z = 1 for the GO and 
Gl simulations. At z = the characteristic mass at the 
high-mass end is roughly similar for the two resolutions, 
around 2.10^^ Mq. The a slope is quite shallow, around 
or higher than —1, whereas the (3 slope is less than —1.5. 
These values are consistent with the ones obtained in dif- 
ferent surveys (Table [2J). It is quite remarkable that such 
a simple galaxy formation model reproduces the shape of 
the galaxy luminosity function. The a slope does not show 
strong evolution between z = and z = 1 , contrary to the 
characteristic mass decreasing as the redshift increases. 
The transitional mass and the normalization parameter 
show clearly the infiuence of resolution: the a power-law 
extends on a higher mass range for the highest resolution, 
changing the transitional mass: Mt ~ 6.10® M© for the 
Gl simulation whereas Mt ~ 2.10^ M© in GO. The bot- 
tom panel in Fig. IA.2I shows a comparison between both 
resolutions at z = 0. Indeed a higher resolution allows 
the formation of lower mass galaxy-like objects. The mass 
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GO 
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-21.8 


-16.6 B band 


Lovedav et al. fl992l 


0, 


.014 


-0, 


.97 


-20.27 




GO 


1.7 


-24.6 


-19.4 K band 


Marzke et al. (19M) 


0, 


.04 


-1, 


.0 


-20.02 












Ellis et al. fl996) 














Gl 


7.4 


-21.1 


-14.5 B band 


T.oveda,v n997^ 


0, 


.0154 


-0, 
-1, 


.94 

.82(/3) 


-20.42 
-14.84(A4t) 




Gl 


3.7 


-23.8 


-17.3 K band 


Zucca et al. ( 1997~l 


0, 


.02 


-1, 


.22 


-20.38 


Table 4. 


Mass-luminosity ratios and characteristic mag- 
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RatcUffe et al. ( 1998) 


0, 


.017 
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.04 


-20.45 












Norbere et al. ('2002) 


0, 


.0168 


-1, 


.21 


-20.43 
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Glazebrook et al. fl995~) 


0, 


.029 


-1, 


.04 


-23.5 




GO (z = 1) 




-21.4 


-16.8 B band 


Gardner et al. f 19971 


0, 


.0166 


-0, 


.91 


-23.9 




GO \z = 1) 




-24.2 


-19.6 K band 


Lovedav (20001 


0, 


.012 


-1, 


.16 


-24.35 










Kochanek et al. f20011 


0, 


.0116 


-1, 


.09 


-24.16 




Gl {z = 1) 




-20.3 


-13.9 B band 


Cole et al. ('20011 


0, 


,0116 


-0, 


.93 


-24.13 




Gl \z = 1) 




-23.0 


-16.7 K band 



Table 2. Fit parameters of the galaxy luminosity func- 
tion extracted from the literature for different surveys us- 
ing a standard Schechter function (Eq. (|lU|l l. Magnitudes 
are computed for h = 0.7. The normalization parameter is 
expressed in fe ^Mpc"^. A two power-law Schechter func- 
tion is used in iLovedaV {imf) wi th a corr ection for the 
/3 parameter in Loy edLay (19981. Zucca et a l. (1997) de- 
scribes the faint part of the galaxy luminosity function by 
a power-law introducing the /? slope and the magnitude 
Mc (line noted (6)). 



GO 

2 = 0.00720 2.53 x 10^ 
2 = 1 0.00824 1.67 X 10^ 



-0.86 2.13 X 10"* 
-0.87 2.55 X 10' 
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Gl 



2 = 1 



0.0146 



1.29 X 10" 



-0.96 
-1.08 



6.35 X 10** 
3.69 X 10* 



-1.65 
-1.59 



-1.52 
-2.02 



Table 3. Fit parameters (Eq. (|ll|l l computed for the 
galaxy-like object mass function (Fig. [SJ for the GO and 
Gl simulations at z = and z — \. The characteristic 
masses are expressed in M© and the normalization pa- 
rameter in /i^Mpc"'^. 



functions tend to be steeper in the intermediate range at 
z = 1, and the characteristic mass decreases. We will 
return to these differences at the end of this section. 

The integration of the galaxy mass function gives 
the mean mass density, jm = / 0(A/)MdM. Using the 
fits displayed in Table |31 we compute jm over the range 
M > 10^ M0. At 2; = this quantity is 8.73 x 10^ and 
1.84 X 10^ h Mq Mpc"^ for the GO and Gl simulations, 
respectively. At z = 1 the mean mass densities are lower: 
6.9 X 10^ and 9.9 x 10^ h Mq Mpc"^ for the same sim- 



Table 5. Characteristic magnitudes in the B and K 
bands for the GO and Gl simulations at z = 1 estimated 
with the mass-to- light ratios at z = given in Table 01 



ulations, respectively. Normalized to the critical density 
parameter at z = 0, pc,o = 2.67 x 10" Mq Mpc"^, the 
stellar density parameter decreases from fi* = 0.00451 to 
0.00357 between z = and z = 1 for the GO simulation, 
and from f]* = 0.00954 to .00511 for Gl. The v alues 
at z = can be compared to iFukueita et al.1 ()l998(l . re- 
porting a central value of = 0.0035 (see their Table 3) 
which accounts for stars in spheroids, disks and irregulars. 

In order to compare the characteristic masses and 
Mt with the characteristic magnitudes of the galaxy lumi- 
nosity function, we use the mean mass densities to derive 
the mass-to-light ratios of the galaxy-like object catalogs. 
We express the mass-to-light ratio as 7 = Jm/jl where jl 
is the mean luminosity density, jl — /o°° (j){L)LdL. Note 
that this ratio is a "stellar" mass-to-light ratio since the 
galaxy-like objects enclose only "stellar" material. The 
mean luminosity density is taken to be 2.5 x lO'^ and 
5.10^ h L(T) Mpc ~'^ in the bj and K bands, respectively 
(ICole et al.ll200l . Table H displays the mass-to-light ra- 
tios in each band and for the two resolutions. Values are 
higher in the B band and also for the highest resolution 
simulation. They are between 1.7 and 7.4, c onsistent with 
observational results (* Fukugita et al.l ()l998ll give mass-to- 
light ratios of 6.5 for spheroids and 1.5 for disks). Recall 
that these mass-to-light ratios are computed over an ex- 
tensive range of mass, and that this quantity is likely to 
depend on mass. 

From these mass-to-light ratios we estimate the 
characteristic magnitudes M* and Mt corresponding to 
the characteristic masses and Mt (in the expression 
M — Mq = —2.5 x log(L/L0) solar magnitude values 
are taken to be Mf^ = 5.3 and M% — 3.3). Results are 
given in Table 01 In the B band M* is brighter than 
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Fig. 9. Dark matter halo mass function at different red- 
shifts for the GO simulation. 



observational values (between -20 and -20.5, see Table[21). 
In the K band is in agreement with observations 
(around -24) . Indeed the infrared band is a better tracer of 
the stellar content, and is not dominated by young stellar 
population strongly emitting in ultraviolet and affected 
by dust extinction. Converting mass into luminosity using 
mass-to-light ratios is then expected to be more reliable 
in the K band. The tran sitional magnitu des defined 
in iLovedavl l(l997l) and Zucca et all l|l997t) are -14.8 
and —17.7 (Table EJ, respectively, and our transitional 
magnitudes in Table 01 are consistent with these values. 
Making the crude assumption that the mass-to-light 
ratio is constant between z = and z — 1, we derive 
the characteristic magnitudes at 2; = 1 corresponding 
to the characteristic masses quoted in Table O Table |S1 
shows that is roughly half a magnitude fainter in 
the B band than at z = and this trend is similar to 
the evolution wi th redshift of in the survey Autofib 
llEllis et al.lll'991 . 

Before going further into the description of the galaxy 
population we compare the galaxy mass function with the 
dark matter halo mass function. Figure shows that the 
mass range extends from 10° to more than lO^'' M©. At 
any redshift the number density of dark matter halos per 
bin of mass increases as the mass decreases. No strong 
decrease at the high-mass end is seen, as in the galaxy 
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Fig. 10. Analytical fits (solid curve) of the dark matter 
halo mass function (cross) for the GO simulation at 2; = 
and z — I, computed using Eq. (|12|) . Dashed curve (iden- 
tical to the solid curve at M > 10^" -^^o) is a three power- 
law analytical fit. 



mass function, at low redshift because of our small com- 
putational v olume, this behav ior being expected at much 
higher mass ij.Ienkins et al.l200lt> . A slight decrease is seen 
around 10^^ Mq and we thus adopt a two power-law fit in 
the range M > 10° M0, such as: 



dN 
dM 



Mt- 



M 



M 



(12) 



where 0* the normalization parameter, a and 7 the slopes 
in the high and low-mass range, respectively, and Mt^ a 
transitional mass between the two power-laws. Fits and 
their parameters at z = 1 and z = are described in 
Fig. Uni (sohd curve) and Table El The Mt, value is 3.4 x 
10" Mq at z = and around 9.10^^ at z = 1. The halo 
mass functions are much steeper than the galaxy mass 
functions in the intermediate mass range. 

However the halo mass functions steepens at the low- 
mass end, M < 10^" Mq, and the analytical function in 
Eq. H12|l does not correctly fit the mass functions for the 
entire mass range. We then use a three power-law func- 
tion by multiplying dN/dM by the term [1 -|- (M/Mt)^], 
similarly to the galaxy mass function. The fit is shown in 
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2 = 


0.000426 


3.40 X lO" 


-1.60 


1.39 








0.000743 


2.95 X 10^^ 


-1.50 


1.40 


7.13 X 10^ 


-1.86 


Z=l 


0.001 


9.02 X 10^^ 


-1.62 


1.57 








0.001 


9.02 X 10^^ 


-1.60 


1.57 


6.11 X 10^ 


-1.89 



Table 6. Fit parameters (Eq. 112|) ') computed for the dark matter halo mass function fFig. llO|l for the GO simulation 
at z = and z = 1. The characteristic masses are expressed in Mq and the normalization parameter in /I'^Mpc^'^. 
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Fig. 11. Conditional mean of the galaxy-like object epoch 
of formation given a mass M , (i|M), for the GO simulation 
catalog at z = (today is at the upper end of the ordinate 
axis). The dispersion ±cr around the conditional mean is 
shown by the error bars. 



Fig.Uniby the dashed curve and Table El gives this second 
set of now 6 parameters. 
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Fig. 12. Top panel : Projected distribution of the galaxy- 
like objects inside a radius < R — 400 h^^ kpc around 
the mass center of three dark matter halos with a mass 
higher than 10^^ M©, for the GO simulation. The halo mass 
centers are at the origin of each plots, scaled in h^^ Mpc. 
Each symbol represents the mass range of galaxies: 10* < 
M < 10^ Mq (cross), 10^ < M < 10" Mq (triangle), 
M > 10^^ Mq (circle). Lower panel: Same as the top 
panel, but for R = 600 h^^ kpc. 



3.3. Epoch of formation 

The galaxy mass function characterizes the galaxy popu- 
lation at a given epoch and does not give any information 
about the background of objects. Catalogs are likely to 
mix galaxies with different properties. The epoch of for- 
mation is then a first insight into their history. Since each 
stellar particle carries a formation epoch a* and a mass 
m*, the epoch of formation of any object is determined 
from each formation epoch of its stellar particles weighted 
by their mass, ^(a,m,)/^m,. 

Figure El displays the conditional mean of the for- 
mation epoch of galaxies given a mass A/, {t\M) (the 
age is the Hubble time minus the epoch of formation), 
computed for the catalog at z = 0. The dispersion around 
the mean is shown by the error bars. This quantity is 
only shown for the GO simulation as the conclusions are 
the same for Gl. The epoch of formation decreases with 
increasing mass, for M > 10^^ Mq galaxies, whereas 
the low-mass galaxies have an epoch of formation glob- 
ally constant. High-mass objects show early epoch of 
formations (see also iPearce et al.l l|200ll) 'l but low-mass 



objects have formed more recently. The low dispersion 
around the conditional mean for the high-mass galaxies 
reflect the fact that their star formation has considerably 
slowed down at low redshift, contrary to the low-mass 
galaxies showing a large dispersion. This suggests that a 
part of the low-mass objects underwent very recently, or 
currently have, a star formation activity at low redshift. 
Another part of these objects have formed, like the 
high-mass objects, at high redshift and have stopped 
their star formation. Our results show that most of the 
stellar mass have formed by redshift z = 1, 68% in the GO 
simulation (37% by redshift z = 2 and 16% by redshift 
z = 3). In the Gl simulation 61% have formed by redshift 
z = 1 and 12% by redshift z = 3. This trend and the 
fact that the most massive systems have quite old epochs 
of for mation are similar a part of ISprineel &: Hernauisll 
lj2003}) conclusions, although they have conducted a 
sophisticated study of the evolution of the star formation 
rate, using a large set of high-resolution simulations 
based on SPH methods and including star formation, 
supernova feedback and galactic outflows. Therefore the 
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galaxy catalog seems to be populated on the one hand 
by high-mass galaxies being in majority early-formed 
galaxies (see Fig. (Sjand Fig. On the other hand, the 
low-mass end of the galaxy mass function is dominated 
either by early-formed galaxies or late-formed galaxies 
with a star formation activity. Separating these two last 
sub-populations would change the shape of the galaxy 
mass function. Regarding the fact that the faint end of 
the galaxy luminosity fun ction is likely to b e populated 
by st ar-forming galaxies ( Zucca et "all Il997t iMadgwic^ 
120021) . this numerical result is quite encouraging. 

The population of high-mass, early-formed galaxies 
could be the observational counterpart of red, passive el- 
liptical galaxies in the center of galaxy clusters. To ad- 
dress this issue we focus on the galaxy-like object distri- 
bution inside the highest mass dark matter halos, with 
M > 10^^ M0. Figure El displays the galaxy distribution 
at a distance < R around the mass center of three dark 
matter halos, randomly chosen in the computational vol- 
ume among the most massive ones. Galaxies are plotted, 
with symbols according to their mass range, inside a ra- 
dius of R — 400 (top panel) or 600 h'^ kpc (lower panel) 
around each mass center halos. Dark matter halos host a 
whole population of galaxies, from M = 10* to more than 
10" Mq. 

It is remarkable that each high-mass halo includes in 
its center a high-mass galaxy-like object (whose position 
at center is not defined a priori). Moreover the biggest 
halo contains the biggest galaxy-like object of the catalog. 
It has been pointed out that the most massive galaxy-like 
objects are also old objects and this result is consistent 
with the observational evidence that galaxy clusters have 
cD type galaxy in their center. It is interesting to note 
that inside a radius of R — 400 h^^ kpc no galaxy- like 
object with intermediate mass, between 10^-10^^ Mq, is 
found in the proximity of the high-mass galaxy, although 
low-mass objects, M < 10^ M©, are present. Galaxies 
of intermediate mass appear when the radius around the 
mass center increases (lower panel). 

4. Photoionization processes 

We now turn to the GP simulation including an ultravi- 
olet background radiation (see Section 12) to analyze the 
influence of photoionization processes on the galaxy-like 
object properties. The same properties as in the previous 
section are discussed. Figure El overplots the cosmic star 
formation rate density for the GP simulation. Although 
the amplitude is lower than for Gl, the general trend is not 
affected by photoionization processes and the SFR density 
decreases from z ^ 2 to the present. Recall that the GP 
simulation has a smaller computational box length than 
GO and Gl and that non-linear long wavelengths are likely 
to be missing at low red shift. A dip appears around z — 6 
l)Barkana fc Loebll200(]|) suggesting that low-mass stellar 
particles do not form above this redshift in low-mass struc- 
tures where gas is now heated by photoionization pro- 
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Fig. 13. Galaxy-like object mass function for the GP 
(cross) and Gl (square) simulations at z = 0. 
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Fig. 14. Same as Fig. but for the GP simulation. 



cesses. Indeed the gas distribution in the temperature- 
density diagram (middle panel in Fig.jSJ shows that these 
processes are dominant over cooling processes only in low 
and middle density regions. 

Figure 1131 compares the galaxy mass function for the 
GP and Gl simulations at z = 0. It reveals the influence of 
photoionization in low-density regions: the most dramatic 
differences are seen at the low-mass end, M < 10^ Mq. 
The slopes and characteristic masses are similar as for Gl. 
The slight decrease at the high-mass end is likely due to 
the smaller computational volume. The formation epochs 
of galaxies are plotted in Fig. 1141 down to 10® M©, below 
this mass the galaxy-mass function decreases. The plot 
presents much more dispersion at M > 10* Mq than in 
Fig. 111! the conditional mean formation epoch varying be- 
tween 4 and 9 Gyr instead around 6 in the GO simulation. 
Galaxies with a mass lower than 10* Mq form at early 
epochs since the formation of low-mass stellar particle 
stops at low redshift (see Fig. IA.l| . Similarly to previous 
results (see Fig. lll|) the highest mass galaxies seem to form 
at early epochs. Our results suggest that, with our choice 
for the F{z) function, the photoionization processes have 
no dramatic influence on galaxy formation over a mass 
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out galaxy formation. 



range higher than 10 '' Mq (see also iQuinn et al.1 l|l996^ ; 
IWeinberg et al More numerical investigations are 

nevertheless required to allow definitive conclusions. 



Fig. 16. Evolution with redshift of the stellar particle 
mean mass (solid lines) with the dispersion icr around 
the mean (dashed and dotted lines) for the GNE (thick 
lines) and Gl simulations (thin lines). 



5. Non-equipartition processes 

We have investigated in Paper I the influence of additional 
dissipative processes on the cosmological plasma: non- 
equipartition processes between ions, neutrals and elec- 
trons may change the gas thermodynamic properties. The 
astrophysical implications are now discussed and quanti- 
fied. 

Paper I shows that the low-density, outer regions of 
gravitationally bound structures are found to be warmer 
in simulations including non-equipartition processes than 
in simulations in which equipartition is forced. This results 
in a warm gas fraction, illustrated by the gas distribu- 
tion in Fig. [3 (bottom panel). Moreover we show that the 
cooling timescale of the warm plasma is longer than the 
cooling timescale of the same regions in simulations with 
forced equipartition. Figure 1151 gives the baryonic mass 
fractions in different ranges of temperature for the simu- 
lation including non-equipartition processes from Paper I, 
the GNE simulation but without galaxy formation, and 
should be compared with Fig for Gl: the fraction of gas 
with a temperature higher than 9.10'^ K is similar in both 
figures but the gas is not distributed in the same phases: 
at z = 5, for example, 3.1% of the plasma constitutes the 
warm phase and 6.4% the cold phase in the simulation 
with non-equipartition processes, whereas these phases are 
0.2% and 9.1% in Gl, respectively. The decrease in avail- 
able cold gas needed for a galaxy to form is expected to 
affect the galaxy formation process. 

We then analyze the GNE simulation, the same sim- 
ulation as performed in Paper I but now including galaxy 
formation (Table As the influence of non-equipartition 
processes is likely to be dominant at epochs before the end 
of the reionization epoch, we do not include photoioniza- 
tion processes and discuss the results at redshifts lower 
than 10 to mimic what happened before this epoch. 



Figure ^] compares the evolution with redshift of the 
stellar particle mean mass between the GNE and Gl sim- 
ulations. The galaxy formation process starts at a lower 
redshift in the former simulation implying a higher mean 
mass. Such differences are also seen in the stellar parti- 
cle mass functions (bottom panel in Fig. IA.1|I . Figure El 
compares the galaxy-like object mass function in GNE 
and Gl. Differences can be seen at the low and high-mass 
ends: at z = 8 the numerical density of galaxies with a 
mass higher than 5.10*^ Mq is 0.2 /i^Mpc"^ in the Gl 
simulation but only 0.08 /i^Mpc"^ in GNE. This repre- 
sents a 60% decrease in the number of objects in the latter 
simulation. Moreover the mass of the biggest galaxies is 
lower in GNE, W against 5.10^ Mq. At z = 6 differ- 
ences are less dramatic and only appear for the low-mass 
end of the galaxy mass function. As discussed in Paper 
I the influence of non-equipartition processes is dominant 
in shallow potential wells making gravitational compres- 
sion unable to heat cosmological plasma to temperatures 
higher than 10^ K. Added to the fact that potential wells 
become deeper as cosmological evolution proceeds, this 
explains the decrease in the fraction of out of equilibrium 
plasma at lower redshift. 

The effects on galaxy properties due to thermodynamic 
modifications of the inter-galactic medium are clearly 
shown here. A follow-up paper (Courty & Alimi, in prepa- 
ration) using the same simulations as in the present pa- 
per will discuss changes in the clustering properties of the 
galaxy-like objects, giving some insights into the physical 
origin of cosmological bias. 

6. Conclusions 

We examine global galaxy properties and the connection 
between galaxy formation and the thermodynamics of the 
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Fig. 17. Galaxy-like object mass function for the GNE 
(square) and Gl (open square) simulations at z = 8 and 
z = Q. 



cosmological gas in Eulerian hydrodynamical simulations. 
The dominant processes known to play a role in galaxy for- 
mation arc included: gravitation, shock heating and cool- 
ing processes. The galaxy formation model consists of the 
identification in the gas distribution of dense and cold re- 
gions. A part of this gas is turned into a stellar particle 
whose collection provides a catalog of galaxy-like objects. 
In addition, the identification of dense regions in the dark 
matter distribution provides a dark matter halo catalog. 
We estimate a number of properties, the cosmic star for- 
mation rate density, the galaxy- like object/dark matter 
halo mass function and the formation epoch of galaxies. 

The most striking result is that considering in the sim- 
ulations these dominant processes gives galaxy properties 
consistent with observations. The cosmic star formation 
rate density shows a peak around z ~ 3 and reproduces at 
low redshift the evolution of the observational star forma- 
tion rate density, namely the sharp decline between z = Q 
and z = 1. More than 60% of the stellar mass has formed 
by redshift z = 1. The galaxy- like object mass function 
presents a significant population of low- mass galaxies and 
shows an evolution with redshift. Moreover, the galaxy 
mass function is well described at low redshift in the range 
M > 10* Mq by an analytical function combining a stan- 
dard Schechter function with a /3 power-law at the low- 



mass end. The fit, parameterized by two characteristic 
masses and two slopes, appears to be in good agreement 
with the observed galaxy luminosity function. The dark 
matter halo mass function is found to be well fitted by a 
three power-law function in the range 10^ < M < 10^^ 
and is steeper than the galaxy mass function in the in- 
termediate mass range. The estimate of the galaxy for- 
mation epoch shows that high-mass galaxies form, in the 
majority, at early epochs. Moreover the galaxy distribu- 
tion around the mass center of the highest mass dark mat- 
ter halos shows that these halos include in their center a 
high-mass, old galaxy-like object. On the other hand the 
low-mass galaxies present a large dispersion around their 
conditional mean formation epoch, suggesting that they 
include recently formed stellar material. 

The galaxy formation model is simple enough for 
the galaxy formation to depend on the thermodynamic 
properties of the baryonic matter and its distribution. 
Introducing the photoionization processes has a dominant 
effect in the low-density regions. This turns into a decrease 
in the density of the lowest mass objects, M < 10^ M0, 
but no drastic change is seen in the galaxy mass func- 
tion at the high-mass end, neither in the galaxy forma- 
tion epoch at intermediate and high mass end, although 
this quantity shows more dispersion. However these results 
could depend on the adopted intensity of the ultraviolet 
background radiation. On the other hand the introduction 
of the non-equipartition processes between the electrons, 
ions and neutrals of the cosmological plasma results in a 
warmer plasma at high redshift and in not too dense re- 
gions, than in simulations in which equipartition between 
species is forced. Hence the longer cooling timescale delays 
star formation. The galaxy-like object mass function then 
shows a decrease in the density of objects at the low-mass 
end. 

Wc also compare two simulations with different com- 
putational box lengths and show that galaxy properties 
depend on the resolution. The amplitudes of the star for- 
mation rate density and the galaxy mass function arc 
higher for the simulation with the highest resolution than 
with the lowest. Indeed adopting a different resolution af- 
fects the history of the accretion and gas cooling inside 
the dark matter potential wells. This results in a change 
in the fraction of available cold gas at a given time, modi- 
fying therefore the galaxy-like object population. We have 
chosen middle-size box lengths and standard resolutions 
to allow the formation of high as well as low-mass objects. 
This work should then be seen as a qualitative analysis of 
the galaxy formation process and no calibration on obser- 
vations at 2: = is adopted. Altogether the galaxy prop- 
erties draw a consistent picture of the galaxy formation 
process and some common features are independent of the 
resolution: the evolution of the cosmic star formation rate 
density, the shape of the galaxy-like object mass function, 
the facts that high-mass galaxies form at early epochs, 
residing in the highest mass dark matter halos, and that 
low-mass galaxies, for some of them, form at later epochs, 
thus showing a star formation activity. 
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Appendix A: The stellar particle mass function 

Figure rKA\ presents the evolution with redshift of the stel- 
lar particle mass function in the three kinds of simulations. 
In Gl the stellar particle mass ranges between 3.10^ and 
2.10^ M0 at z = 0, typically that of the globular clusters. 
We note that the mass range of the particles increases as 
the redshift decreases: at any epoch, not only have low- 
mass stellar particles formed, but also larger ones. 

In the GP simulation (middle panel) the formation of 
the low-mass stellar particles is stopped when photoion- 
ization processes become dominant, for redshifts less than 
7 (after the epoch of the sharp steepness of the ultraviolet 
background radiation spectrum, Fig.^. At lower redshift 
the mass function presents an evolution for particles with 
mass higher than 3.10^ Mq and the general shape seen in 
the top panel is modified. 

The bottom panel displays the mass functions for the 
simulation with non-equipartition processes: at high red- 
shift the general shape of upside down "V" seen for the Gl 
simulation, with a similar maximum (around 2.10^ Mq in 
Gl) shows that the numerical density of stellar particles 
is much lower in the GNE simulation than in Gl. 

Appendix B: Sensitivity to parameters 

We briefly discuss how the galaxy mass function is sensi- 
tive to the parameters involved in determining the galaxy- 
like objects and to the resolution. 

The top panel of Fig. IA.2I shows the influence of the 
link parameter 77 used to group stellar particles and ini- 
tially flxed to 0.2 (see Section O. Taking 77 = 0.1 results 
in a steeper (3 slope. The catalog then includes a larger 
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Fig. A.l. Evolution with redshift of the stellar particle 
mass function for the Gl, GP and GNE simulations (from 
top to bottom panels, respectively). The mass function is 
shown at z =9, 8, 7, 6, 4, 2, for the Gl and GNE 
simulations and at z =8, 7, 5, 4, 2, for the GP simulation 
(alternatively filled and open circles, from bottom to top). 



number of low-mass objects. On the other hand, choosing 
rj = 0.3 would result in a shallower (3 slope. Nevertheless 
there is no modification in the a power-law range. 

The middle panel of Fig. IA.2I compares the galaxy 
mass function at z = in the GO simulation consider- 
ing A/min — 1 and Afmin — 10. The parameter Mmin = 1 
means that even a single stellar particle is identified as 
a galaxy-like object. This extends the mass function to- 
wards lower mass, shifting the decrease at the low mass 
end, but does not change the rest of the mass function. 
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Fig. A. 2. Galaxy-like object mass function for the GO 
(-^box = 32) and Gl (ibox = 16) simulations at z = 
computing by considering different values of rj (top panel) 
and Minin (middle panel). The bottom panel compares 
both resolutions. 



Finally the bottom panel of Fig. IA.2I illustrates the 
influence of the resolution. As already seen with the fits 
in Fig. IHl the a and f} slopes and the characteristic mass 

are roughly similar for the two simulations. But in the 
high-resolution simulation, Gl, the transitional mass Mt 
and the decrease at the low mass end are shifted towards 
lower mass, this simulation being able to form lower mass 
objects. 



